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STUDIES  OF  ARTIFICIAL  CRACKS  AT  ROCKET  MOTOR 
LINERS  BY  THE  FROZEN  STRESS  METHOD 
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ABSTRACT 

A  preliminary  investigation  limited  to  stress  freezing  materials  was  undertaken  in  order  to  develop  a 
test  method  for  estimating  the  stress  intensity  factor  (SIF)  distribution  along  the  border  of  cracks  in  the 
bond  line  of  a  bi-material  specimen.  It  was  desired  to  measure  the  influence  of  residual  stress  and  mixed 
mode  effects  for  such  specimens.  A  single  edge  notched  specimen  was  selected,  and  a  test  procedure  was 
devised  for  isolating  the  residual  stress  effect  at  the  bond  line  for  Mode  I.  Subsequently,  mbced  mode  tests 
were  conducted  in  the  absence  of  externally  induced  residual  stress  and  the  SIF  distributions  for  Ki  and 
K2  were  obtained.  It  was  found  that  the  level  of  the  SIF  values  were  increased  by  the  bond  line  in  cracked 
homogeneous  bondline  specimens  for  bondlines  nearly  normal  to  the  applied  load  but  that  this  effect  was 
not  increased  in  bonded  bimaterial  specimens.  For  larger  angles  of  inclination  to  a  normal  to  the  applied 
load  no  increase  in  the  SIF  above  the  homogeneous  result  was  found. 

INTRODUCTION 

When  a  crack  is  located  at  the  interface  between  two  materials,  near  tip  stress  analysis  is  complicated 
by  the  fact  that  the  mode  1  and  mode  2  stress  intensity  factors  (SIF)  cannot  be  decoupled  [1].  However,  if 
the  materials  are  incompressible,  and  plane  strain  exists,  the  near  tip  interface  equations  reduce  to  the  more 
familiar  form  for  isotropic  and  homogeneous  elastic  solids  [2]. 

When  cracks  occur  between  a  solid  rocket  motor  and  its  rubber  liner,  the  above  conditions  are  approx¬ 
imately  realized  and  occur  at  a  rather  low  modulus  ratio  between  the  two  materials.  The  first  author  and 
his  colleagues  have  had  considerable  success  in  measuring  three  dimensional  effects  in  cracked  isotropic  and 
homogeneoiis  bodies  using  a  refined  frozen  stress  method  [3].  The  frozen  stress  method  is  briefly  summarized 
in  Appendix  I.  In  order  to  explore  the  feasibility  of  such  an  approach  in  measimng  three  dimensional  and 
residual  stress  effects  in  such  cracks  as  described  above,  a  series  of  experiments  were  conducted  on  thick 
plates  in  tension  containing  bondline  cracks  between  two  different  incompressible  materials. 

The  experiments  can  be  divided  into  two  parts: 

i.  A  study  of  externally  induced  residual  stress  effects  on  a  bond  line  crack  under  Mode  I 
loadings  and 

ii.  A  study  of  mixed  mode  effects  in  the  absence  of  externally  induced  residual  stresses. 

SPECIMEN  CONFIGURATION.  MATERIALS  AND  TEST  PROGRAM 
The  specimen  configuration  used  in  all  tests  is  pictured  in  Fig.  1.  All  loads  were  parallel  to  the  long 
dimension  of  the  specimens.  The  materials  and  their  properties  au:e  given  in  Table  I. 


Table  I 


Material 

^critical 

Enot 

fr 

Araldite 

240'’F 

2698  psi  (18.60  MPa) 

1.64  psi-in  (286.9  Pa-m) 

Aral- Alum 

240°F 

5349  psi  (36.88  MPa) 
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The  thermal  coefficient  of  expansion  for  both  materials  was  matched  exactly  at  15.3  x  10  ®  per  at 
room  temperature  and  119  x  10“®  ±  20  xlO""®  per  at  240®F.  The,  Araldite-aluminum  specimens  were 
produced  in  opaque  form  by  adding  25%  by  weight  of  aluminum  powder  to  the  pure  transparent  araldite, 
a  common  photoelastic  material.  The  bonchng  agent  was  a  commercial  glue  PLM-9  made  by  Photoelastic, 
Inc.  Raleigh,  NC.  The  test  specimen  materials  were  supplied  by  Survey  Technology  Centre,  London. 

The  testing  program  involved  three  test  specimens  for  each  crack  orientation:  one  homogeneous  specimen 
(B-1),  one  bonded  homogeneous  specimen  (B-2)  and  one  bonded  bimaterial  specimen  (B).  This  was  done  in 
order  to  isolate  residual  bond  line  effects  from  bimaterial  effects. 

ANALYTICAL  PROCEDURES 

In  order  to  provide  an  analytical  control,  the  following  procedure  was  used  to  estimate  theoretical  values 

of  and  and  ^ 

The  values  of  Ki^th)  and  K2{th)  were  based  on  two  dimensional  solutions  for  edge  cracks  m  semi-mnmte 
plates  [4],  [5]  adjusted  to  finite  width  plates  using  [6],  [7]  and  correcting  for  the  notch  effect  using  [8].  In 
adjusting  the  infinite  plate  solution  to  the  finite  plate  solution  for  inclined  cracks  #  0),  the  projection  of 
the  crack  normal  to  the  direction  of  the  applied  load  was  used. 

Experimental  results  were  obtained  from  frozen  stress  photoelastic  experiments  and  optical  data  were 
converted  into  stress  intensity  values  using  the  algorithms  and  procedures  described  in  Appendix  H. 

TEST  RESULTS 


RESIDUAL  STRESS  TESTS 

In  order  to  generate  a  residual  stress  field,  thr-ee  specimens  B-1,  B-2  and  B  with  /?  =  0  (Fig.  1)  were  cycled 
to  280°F  to  exacerbate  the  thermal  coefficient  variations  at  elevated  temperature  between  the  araldite  and 
the  araldite-aluminum  materials.  Then  after  stress  freezing  using  the  280®F  maximum  temperature,  thin 
slices  were  removed  mutually  orthogonal  to  the  crack  front  and  the  crack  plane  (L,M,R  in  Fig.  1)  and  the 
slices  were  analyzed  photoelastically  (Ap  I)  and  SIF  values  were  determined  for  each  slice  (Ap  H).  Fringe 
patterns  from  the  cracked,  homogeneous  specimen  were  as  expected  and  are  not  shown  here.  Fig.  2  shows 
the  no  load  (after  thermal  cycling),  after  stress  freezing  (through  thickness)  and  5th  multiple  of  a  side  slice 
for  the  B-2  specimen  and  Fig.  3  shows  similar  results  for  the  B  specimen.  Clearly,  the  B  specimen  shows 
higher  no  load  residual  stresses  and,  in  all  cases,  a  fringe  pattern  was  produced  near  to  the  bond  line.  As 
will  be  seen  later,  the  increase  in  the  Ki  due  to  the  bi-material,  or  modulus  effect  is  believed  due  to  the 
residual  stresses  produced  by  thermal  cycling. 

TESTS  WITHOUT  INDUCED  RESIDUAL  STRESS 

Frozen  stress  tests  without  the  thermal  cycling  to  induce  residual  stresses  were  conducted  for  values  of  (3 
of  45°,  30°,  15°  and  0:  Fringe  patterns  from  the  homogeneous  bonded  tests  of  0°  and  30°,  as  shown  in  Fig. 
4.  For  =  30°  in  the  through  thickness,  no  load  case,  two  fringes  are  visible  along  the  bond  line  (Fig.  4a). 
Machining  stresses  along  notch  sides  tend  to  anneal  out  during  stress  freezing.  Upon  loading,  some  twenty 
fringes  appear  locally  outside  the  bond  line  (Fig.  4b).  Fig.  4(c)  shows  the  pattern  from  a  slice  less  than 
1mm  thick.  The  dark  lines  parallel  to  the  bond  line  denote  bond  line  araldite  interfaces. 

Fringe  patterns  for  the  0°  case  are  shown  in  Figs.  4d  through  4f.  The  no  load  and  loaded  through 
thickness  patterns  are  similar  to  the  30°  results  qualitatively  as  is  the  slice  pattern.  The  no  load  fringes 
parallel  to  the  bond  line  result  from  bond  curing. 

In  the  15°  homogeneous  bonded  test,  a  crack  emanated  from  the  notch  tip  under  load  at  the  notch  tip, 
ran  to  a  bond-araldite  interface  and  proceeded  to  grow  along  the  interface. 

The  normalized  SIF  results  from  all  tests  are  tabulated  in  Table  II.  Several  general  comments  can  be 
made. 

i)  All  Fj'h  values  are  two  dimensional  analytical  results.  {Fij)th  values  axe  estimated  from  {Fi)th  analysis. 

ii)  The  F/  distributions  through  the  thickness  revealed  a  slight  elevation  in  the  SIF  near  mid-thickness  as 
expected  in  every  case  except  the  45°B-M  result. 

hi)  Only  the  15°  and  the  0°  cases  revealed  a  significant  increase  in  F\I{Fi)tHavg  bimaterial  specimens 

(14%,  16%).  However,  this  same  order  of  increase  was  found  for  the  0°  homogeneous  bonded  specimen 
(19%)  and  would  be  expected  for  the  15°  case  as  well. 

iv)  Values  of  {Fn I{Fu)th)avg  were  generally  higher  than  their  {Fj I{Fi)th)  avg  counterparts  but,  since 
we  do  not  have  an  accurate  two  dimensional  analytical  result  for  this  case,  no  particular  significance  is 
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Table  H 


=  2D  Theoretical  estimates 

=  Kj'd^fm  ;  F^  =  K^l'd>!m 

a  =  0.375  in.  (9.53  m.m.)  [Fig.  1] 

B1  =  Homogeneous  specimen 

B2  =  Homogeneous  bonded  specimen 

B  =  Bi-material  specimen 

LR  =  Average  of  left  edge  and  right  ecfee  slices  [Fig.  1] 

LMR  =  Average  of  left  edge,  middle,  and  tight  edge  slices  [Fig.  1] 
For  projectioiis  of  mode  1  and  mode  2  results  are  used. 


Specimen 


Bl-LR 

B1-M,LMR 

B2-LR 

B2-M,LMR 

B-LR 

B-M,LMR 


Bl-LR 

B1-M,LMR 

B2-LR 

B2-M,LMR 

B-LR 

B-M,LMR 


Bl-LR 

B1-M,LMR 

B-LR 

b-m,lmr 


Bl-LR 
Bl-M,  LMR 
B2-LR 
B2-M,  LMR 
B-LR 

B-M.LMR 


F, _ ^ _ 

/?  =  0°  (with  no-load  thermal  cycle) ;  [Fig.  1] ;  (f]  );„  =  1 .5 1 

_ 5=  =  7.00psi.  (0,048  Mpa.) _ 

1.70  113 

1.64,  1.68  -  109,  1.11 

1.78  118 

1.71,  1.75  -  1.13,  1.16 

1.81  -  1.20 

1.91,  1.84 _ _  _ 1.26,  1.22 _ 

(without  no-load  thennal  cycle) ;  [Fig.  1] ;  (^  )7h  ~  1-^1 

_ _ g  =  8.75  psi.  (0.060  Mpa.) _ _ 

1 37  0.907 

1.67,1.47  -  1.11,0.974 

1.73  1.15 

1.95,  1.80  -  1.29,  1.19 

1.72  1.14 

1.82,  1.75 _ _ 1.21,  1.16 _ 

^  =  15';  [Fig.  1] ;  (i^)^  =  1.41 ;  =0.181 

=  14.0  psi.  (0.097  N^a.) ;  =  8.75  psi.  (0.060  MPa.) 


1.31 

1.62,  1.42 
1.52 

1.79,  1.61 


0.179 

0.221, 0.193 
0.209 

0.244,0.231 


0.929 
1.15,  1.01 
1.08 

1.27,  1.14 


y0  =  3O”;  [Fig.  1]  ■AF^)r„  =  =0.339 

=  7.00  psi.  (0.048  Iv^a.) ;  Ob2,b  -  ^*^  0 1”^-  (0  097  MPa.) 

1.03  0.320  0.888 

1.14,1.07  0.353,0.330  0.983,0.922 

1.05  0.323  0.905 

1.16.1.08  0.357,0.334  1.00,0.931 

1.02  0.317  0.879 

1.16,1.07  0.357,0.330 _ 1.00,  0.922 _ 

^  =  45”;  [Fig.  1):  (/=;)„,  =  0.837  ;  (F,  )j„  =0.432 
g  =  7.00  psi.  (0.048  Mpa.) 


Fj{F,)r 


0.989 
1.22,  1.07 
1.15 

1.35,  1.28 


0.944 
1.04, 0.973 
0.953 

1.05,  0.985 
0.935 
1.05.  0.973 


Bl-LR 
Bl-M.  LMR 
B2-LR 
B2-M.  LMR 
B-LR 

B-M.  LMR 


0.856 

0.921,0.877 

0.843 

0.873.  0.852 
0.893 

0.857.  0.881 


0.494 

0.532,  0.506 
0.486 

0.504,  0.492 
0.516 

0.494.  0.508 


1.02 

1.10,  1.05 
1.01 

1.04,  1.02 
1.07 

1.02.  1.05 


1.14 

1.23.  1.17 
1.13 

1.17,  1.14 
1.19 

1.14,  1.18 
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attached  to  these  elevations. 

v)  The  residual  stress  0®  results  from  the  bimaterial  specimen  were  similar  to  the  0®  result  for  the  bimaterial 
test  without  thermal  cycling  suggesting  that  the  mismatch  in  thermal  coefficients  between  the  araldite 
and  the  Al-araldite  above  critical  temperature  has  little  effect  on  the  resulting  K\  determination. 

SUMMARY 

A  series  of  tests  were  conducted  in  order  to  provide  a  data  base  for  the  development  of  a  test  for  evaluating 
stress  intensity  factor  distributions  along  bond  lines  between  incompressible  materials  of  different  Moduli  of 
Elasticity.  It  was  desired  to  separate  the  bond  line  effect  from  the  bi-material  effect. 

In  a  tot  test  series,  residual  stress  was  induced  prior  to  live  loads,  and  fringe  patterns  from  homogeneous, 
homogeneous  bonded,  and  bimaterial  specimens  were  examined.  Table  m  summarizes  the  results. 


Tkble  III 


Effect  on  (from  slicine)  of 

B-1 

Machining  Stresses  (MS) 

11 

B-2 

Machining  Stress  +  Bond  Line  Stress  (BS) 

16 

B 

Machining  Stress  +  Bond  Line  +  Modulus 

22 

As  shown  in  Table  H,  however,  variations  in  the  SIF  across  the  plate  thickness  were  small  and  of  the  order 
of  the  accuracy  of  the  method  of  analysis  (3.6%).  However,  they  did  show  the  expected  distribution  in  the 
bimaterial  specimen.  Moreover,  the  procedure  allows  separation  of  bond  stress  and  bimaterial  stress  effects. 

Tests  without  induced  residual  stress  were  conduct^  on  test  specimens  containing  bond  line  artificial 
cracks  inclined  to  the  load  direction  at  the  angles  of  0®,  15®,  30®  and  45®.  The  modulus  ratio  in  the  bimaterial 
specimens  was  1.97.  Stress  intensity  factors  determined  at  intervals  through  the  specimens  thicknesses  by 
the  frozen  stress  method  were  normalized  with  respect  to  two  dimensional  solutions  for  Mode  I  and  two 
dimensional  estimates  for  Mode  n. 

As  a  result  of  these  studies,  it  was  found  that: 

i)  For  values  of  =  30®,  45®  (Fig.  1)  no  increase  in  the  SIF  values  were  observed  due  either  to  the  bond 
line  or  the  bimaterial  effects. 

ii)  For  values  of  =  0®,  15®  (Fig.  1),  an  increase  in  the  SIF  level  occurred,  in  both  cases  of  16%  and  14% 
respectively  in  the  bimaterial  specimens  but  results  indicate  that  all  of  this  increase  is  due  to  the  bondline 
effect,  not  the  modulus  ratio. 

iii)  Mode  n  results  were  larger  than  predicted  by  approximate  analysis  but  again  appeared  due  solely  to  bond 
line  effects. 

iv)  In  all  cases  except  one  the  peak  SIF  occurred  at  mid  thickness.  In  the  single  exception,  ()9  =  45®) 
the  variation  in  the  SIF  across  the  thickness  appeared  to  be  of  the  order  of  the  experimental  scatter, 
(t.e.  <  ±3%). 

v)  It  is  of  interest  to  note  that  there  was  no  effect  on  Ki  of  either  the  homogeneous  bondline  or  the  bimaterial 
specimen  for  values  of  ^  —  30®  and  45®.  This  is  believed  due  to  the  fact  that  the  no  load  bondline  fringes 
were  nearly  parallel  to  the  load  direction  near  the  crack  tip  and  thus  made  no  contribution  to  the  fringe 
gradient  in  the  load  direction  which  is  used  to  determine  K\, 

vi)  For  both  =  0®  cases  and  for  ^  ss  15®,  a  significant  elevation  in  Ki  occurred  in  the  bimaterial  specimens. 
(14%  to  16%).  However,  this  increase  also  occurred  in  the  homogeneous  bonded  specimen  in  the  /?  =  0 
test  without  residual  stress,  also  to  a  lesser  extent  in  the  /?  =  0  homogeneous  bonded  residual  stress  test, 
and  is  expected  in  the  0  =  15®  homogeneous  bonded  test.  The  suggestion  is  that  the  elevation  of  Ki  is 
due  to  the  bondline  stresses  and  not  due  to  the  modulus  ratio  in  the  bimaterial  specimen.  Examination 
of  the  no  load  bondline  fringes  reveals  that  they  are,  indeed  nearly  normal  to  the  load  direction  near  the 
crack  tip. 

In  conclusion,  it  appears  that  using  the  three  specimen  test  procedure,  reasonable  results  can  be  ob¬ 
tained  for  SIF  distributions  through  the  thickness  of  bimaterial  specimens  containing  bond  line  cracks  in 
incompressible  materials.  Tests  employing  this  method  are  being  conducted  on  cracks  parallel  to  the  bond 
line. 
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appendix  I  -  FROZEN  STRESS  PHOTOELASTICITY 
When  a  transparent  model  is  placed  in  a  circularly  polarized  monochromatic  light  field,  and  loaded,  dark 
fringes  will  appear  which  are  proportional  to  the  applied  load.  These  fringes  are  called  stress  fringes  or 
isochromatics  and  the  magnitude  of  the  maximum  inplane  shear  stress  is  a  constant  along  a  given  fringe. 

Some  transparent  materials  exhibit  mechanical  diphase  dmracteristics  above  a  certain  temperature,  called 
the  critical  temperature  (Tc).  The  material,  while  still  perfectly  elastic  will  exhibit  a  fringe  sensitivity  of 
about  twenty  timps  the  value  obtained  at  room  temperature  and  its  modulus  of  elasticity  will  be  reduced 
to  about  one  six  hundredth  of  its  room  temperature  value.  By  raising  the  model  temperature  above  Tc, 
loading,  and  then  cooling  slowly  to  room  temperature,  the  stress  fringes  associated  with  Tc  will  be  retained 
when  the  material  is  returned  to  room  temperature.  Since  the  material  is  so  much  more  sensitive  to  fringe 
generation  above  Tc  than  at  room  temperature,  fringe  recovery  at  room  temperature  fiinge  upon  unloading 
is  negligible.  The  model  may  then  be  sliced  without  disturbing  the  “frozen  in”  fringe  pattern  and  anzJj^ed 
as  a  two  Himpnainnal  model  but  containing  the  three  dimensional  effects.  In  the  use  of  the  method  to  make 
measurements  near  crack  tips,  due  to  the  need  to  reduce  loads  above  critical  temperature  to  preclude  large 
local  deformations,  and  the  use  of  thin  slices,  few  stress  fringes  are  available  by  standard  procedures.  To 
overcome  this  obstacle,  a  refined  polariscope  (Fig.  I-l)  is  employed  to  allow  the  tandem  use  of  the  Post  [9] 
and  Tardy  [10]  methods  to  increase  the  number  of  fringes  available  locally. 

In  fringe  photographs,  integral  fringes  are  dark  in  a  dark  field  and  bright  in  a  bright  field.  Bright  fields 
are  used  throughout  this  paper  for  clarity. 

APPENDIX  II  (ALGORITHMS) 

MODE  T  ALGORITHM  FOR  HOMOGENEOUS  CASE 

Beginning  with  the  Griflath-Irwin  Equations,  we  may  write,  for  Mode  I, 

(Tij  =  (i-j-  =  n,  z)  (11 .1) 

where: 
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aij  are  components  of  stress 
Ki  isSIF 

r,  0  are  measured  from  crack  tip  (Fig.  II-l) 
(T^j  are  non-singular  stress  components. 

Then,  along  0  =  7r/2,  after  tnmcating  aij 


where  r°  =  and  is  constant  over  the  data  range 

Kap  =  apparent  SIF 
Tnz  =  TTiayininm  shear  stress  in  nz  plane 

.  Kap  _  Ki  ^  y/Sr^  /r \  i 
a(7ra)i  a(7ra)i  ^ 


{II-2) 


(J7.3) 


where  (Fig.  11.1)  a  =  crack  length,  and  a  =  remote  normal  stress 
i.e.  vs.  is  linear. 

a{Tra)^  V  ^ 

Since  from  the  Stress-Optic  Law 
{Tnz) max  “  where 

n  =  stress  fringe  order 
/  =  material  fringe  value 
t  =  specimen  (or  slice)  thickness 
and  from  Eq,  11.2 

=  ('^n2)max(87rr)i  =  -5^(87rr)i,  .  .  ^  X'U  *  xu  f 

then  Kap  (through  a  measure  of  n)  and  r  became  the  measured  quantities  from  the  stress  fringe  pattern  at 
different  points  in  the  pattern. 

A  typical  plot  of  normalized  Kap  vs.  y/r^  for  a  homogeneous  specimen  is  shown  in  Fig.  11.2. 

Mixed  Mode  Algorithm 

The  mixed  mode  algorithm  was  developed  (see  Fig.  11-3)  by  requiring  that: 


Hm 

rm->0 


=  0-77-4 


which  leads  to: 


- II -5 

3 


By  measuring  0^  which  is  approximately  in  the  direction  of  the  applied  load,  K2IK1  can  be  determined. 
Then  writing: 

fn_  K\p 

2t  (87rr)2 


max  ^  til  ^ 


one  may  plot  vs. 

o-(7ra)  2 

Knowing,  K*^  K2IK1  and 


\frfa  as  before,  locate  a  linear  zone  and  extrapolate  to  r 
0^,  values  of  Ki  and  K2  may  be  determined  since. 


o  to  obtain  K*. 


K’  =  [(7^1  sin©^  +  2K2COsQlf  +  - 77-6 

Details  are  found  in  Ref.  3.  A  plot  of  vs.  y/rja  will  yield  a  linear  zone  from  which  K'  can 

be  extracted.  Knowing  K*  and  0^,  Ki  &  K2  can  be  determined  from  Eqs.  11-5  and  II-6. 
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c)  5th  Multiple  Side  Slice 


c)  5th  Multiple  Side  Slice 


Fig.  2.  Homogeneous  Bonded  Specimen  (B-1)  with 
Residual  Stress  (Bright  Field)  MF  =  8.75 


Fig.  3.  Bi-Material  Specimen  (B)  With  Residual 
Stress  (Bright  Field)  MF  =  8.75 


a)  No  Load  Thru  Thickness 


d)  No  Load  Thru  Thickness 


c)  Frozen  edge  slice  f)  Frozen  edge  slice 


Fig.  4  Fringe  Patterns  -  homogeneous  bonded  specimen  bright  field  Figs,  a-c,  P  =  30°,  Figs,  d-f,  =  0°  M.F 
=  8.75  for  all  figures. 


